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ABSTRACT 

We investigate whether the subhalos of ACDM galaxy halos have potentials consis- 
tent with the observed properties of Milky Way satellites, particularly those with 
high-quality photometric and kinematic data: Fornax, Leo I, Sculptor, Sextans, and 
Carina. We compare spherical models with isotropic velocity dispersion tensors to the 
observed, circularly averaged star counts, line-of-sight velocity dispersion profiles and 
line-of-sight velocity distributions. We identify subhalos within the six high-resolution 
dark matter halos of the Aquarius Project for which the spherically averaged potentials 
result in excellent fits to each of the five galaxies. In particular, our simple one-integral 
models reproduce the observations in the inner regions, proving that these data are 
fully consistent with ACDM expectations and do not require cored dark matter dis- 
tributions. For four of the five satellites the fits require moderately cusped stellar 
density profiles. The star count data for Leo I, however, do require a cored distribu- 
tion of star counts. Current data suggest that these five satellites may be hosted by 
ACDM subhalos with maximum circular velocities in the range 10 to 30 km s~^. 
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1 INTRODUCTION 

The internal dynamics of the dwarf spheroidal (dSph) satel- 
lite galaxies of the Milky Way (MW) offer perhaps the best 
prospects for investigating the properties of the dark mat- 
ter in the nearby universe. These galaxies are dominated by 
dark matter and the brightest of them are sufficiently close 
that high-precision line-of-sight velocities for large sam- 
ples of their stars can be measured using high -resolution, 
multi-object spectroscopy (| Walker et al.l l2009l \ The cur- 
rent datasets repres ent an improve ment upon the ear liest 
such observations l|Aaronso n 1983 ; iMateo et al.l Il993i ) by 
factors of tens to thousands. Analysis of the kinematic data, 
in combination with improved photometric measurements, 
have not only confirmed earlier indications that the classical 
dSphs are dark matter dominated, but have also revealed the 
surprising property that their mean central densities are sim- 
ilar eve n though their luminosities span a very wide range of 
values jMateo" 19981: 1 Walker et al. l l2007l : lGilmo"re" et al. I l2007l : 
IStrigariet al. 20oi). 

Current cosmogonic theory makes strong predictions for 
the internal structure of dark matter halos. N-body simula- 
tions of halo formation in hierarchical clustering cosmogo- 
nies have shown that halos develop strongly cusped den- 



sity profiles which are almost independent of ha l o mas s 
and cosmological parameters (|Navarro et al ] Ii996l . liggTl i. 
Subsequent simulat i ons have confirme d this result (e.g. 
iNavarro et all I2OI0I : IStadel et all |2009| ). showing, in ad- 
dition, that cuspy profiles are retained even after halos 
fall into larger ones and suffer extensive tidal stripping 
(|Kazantzidis et al.|[200^ : ISpringel et al. 200 8'). Thus, to the 
extent that the dark matter distributions in the inner parts 
of halos have not been significantly disturbed by the galaxies 
forming within them, halo profiles offer a strong and direct 
test of the ACDM cosmogony in a regime not probed by 
microwave background and large-scale structure data. 

An in-depth analysis of the central regions of Milky 
Way satellites is particularly important given that a num- 
ber of recent studies of the structure of these galaxies 
have claimed that shallow central density profiles provide 
a better description of their dark ma tter halos than the 
cuspy profiles characteristic of ACDM ("Goerdt et al.|[2006l : 
ISanchez-Salcedo et al] I2OO6I : ICilmore et al. 200?! ). Models 
with dark matter core radii of ~ 100 pc have been shown 
to provide good fits to the kinemat ic data sets of the classi- 
cal dSphs j Angus fc Diaferidl2009l ). If confirmed, the shal- 
low cores suggested by these studies might indicate a lower 
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central phase-space density than e xpected if the dark mat- 
ter is a cold coUisionless particle (|Tremaine fc Gunnlll979l : 
iHoean fc Dalcanto n"2000\ 

Kinematic studies typically treat dSph galaxies as 
spherical, dynamically equilibrated systems. With these sim- 
ple assumptions there is a strong degeneracy between the 
statistics of stellar orbits (i.e. whether velocity dispersions 
are isotropic, or are radially or tangentially biased) and the 
shape of the stell ar and dark matter density profiles (e.g. 
lEvans et alllixil ). This ambiguity is reflected in the broad 
range of models used in recent attempts to constrain the 
dark matter density profiles of the dSphs ( Strigari et al.l 
l2008l : lLokasll2009l : IWalker et al1l2009l : I Wolf et aiBoogP . The 
intrinsic parameter degeneracies of the models cast doubt 
on the robustness of inferences favoring cored or cuspy cen- 
tral density profiles, even given the high-quality data that 
are now available (Strigari ot al. 2006; Walker ct al. 2009). 
Breaking these degeneracies may only be possible by ex- 
ploiting additional observational constraints, for example, 
through measurement of internal stellar proper motions 
l|Wilkinson et al.ll2002l : IStrigari et aLlbOOTh . 

Motivated by the simple question of whether the ob- 
served dSphs are kinematically consistent with the ACDM 
theory of structure formation, we here use six high- 
resolution ha lo simulations perform ed as part of the Aquar- 
ius Project jSpringel et al.l |2008| ) to search for subhalos 
whose properties would allow them, in principle, to host 
the well observed dSph satellites of the Milky Way. We 
restrict our attention to the five satellites with abundant, 
high-quality stellar kine matic data: Fornax, Sculptor, Leo 
I, Ca rina, and Sextans (|Mateo et al. 1 l2008l : IWalker et al.l 
|2009| ). Assuming isotropic velocity dispersions, we identify 
an Aquarius subhalo with a dark matter potential which 
results in a good simultaneous fit to each satellite's photo- 
metric and kinematic data. 

In addition to focusing specifically on ACDM and mak- 
ing direct use of realistic potentials from the Aquarius sim- 
ulations, our modeling differs from previous work in that we 
simultaneously fit both photometry and kinematics. We al- 
low for mildly cusped steZZar profiles with ~ r~" near the 
centre, where a is in the range to 1. Projections of such pro- 
files are fitted to the photometric data, and with the results 
in hand we make predictions for the kinematic data, both 
traditional second moment (line-of-sight velocity dispersion) 
profiles and full line-of-sight velocity distributions. The lat- 
ter comparison allows us to test whether simple, spherical, 
isotropic, ACDM-based models are consistent with higher 
moments of the observed velocity distribution. 

For each of the five satellites we study, the best-fitting 
Aquarius subhalo provides an excellent statistical fit to the 
data. For four of the five this requires a mildly cusped stel- 
lar distribution similar to those found in brighter early-type 
galaxies. The only galaxy that requires a true core in the 
star distribution is Leo I, but with such a profile its kine- 
matics are still consistent with a ACDM subhalo. We present 
circular velocity curves for the best-fitting subhalo hosts for 
each of the MW satellites. Not surprisingly, for a given satel- 
lite, the circular velocity curves of "good" subhalos are very 
similar at the radii that are well sampled by the stellar trac- 
ers. This is a consequence of our assumptions of spherical 
symmetry and isotropy. Finally, we determine the mass of 
the best-fitting subhalos, both at the time of accretion onto 



the host halo and at high redshift, and we show that these 
quantities have more scatter than the present-day central 
potentials or maximum circular velocities. 



2 THEORETICAL MODELING 

The goal of our theoretical modeling is ultimately to com- 
pare to the full observed line-of-sight velocity distribution for 
each satellite. We begin by discussing how to model velocity 
dispersion profiles and then move on to modeling of the full 
line-of-sight velocity distribution. In our later analysis, we 
will use predictions for line-of-sight velocity dispersion pro- 
files to identify possible subhalo hosts for each satellite, and 
then use the corresponding velocity distributions to check 
that our isotropic models are indeed consistent with the ob- 
served kinematics. 



2.1 Velocity dispersion 

We assume spherical symmetry both for the potential and 
for the distribution of stars within it. The line-of-sight ve- 
locity dispersion of a dSph satellite at projected radius, R, 
can then be written as 



o-;os(-R) 



1 



r>2 



- J?2 



-dr. 



(1) 



where p*(r-) is the stellar density profile and It{R) its two- 
dimensional projection; P = 1 — a^/a'^ is the anisotropy 
parameter, where o-t{r) is the one-dimensional tangential 
velocity dispersion of the stars and (t^('') the corresponding 
radial velocity dispersion. These quantities satisfy the radial 
Jeans equation: 



d(p*(j^ 
dr 
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(2) 



The total mass distribution, M{r), is the sum of the mass 
distributions of dark matter, Mdm{r), and stars, M^{r). 

For our analysis, we use the dark matter mass distri- 
butions, Mdm{r), of subhalos in the next-to-highest resolu- 
tion set of Aqu arius simulations of ga lactic halos (level 2 in 
the notation of ISpringel et al] l|2008l) ). At the present day, 
each of the six halos contains about 45000 resolved subhalos 
within ~ 400 kpc of its centre (corresponding to the radius, 
rso, of the sphere of mean overdensity 50 times the critical 
value) with more than 20 particles, each of mass ~ 10* Mq. 
The structure of each subhalo is well resolved down to a 
physical radius of ~ 100 pc, corresponding approximately 
to tw ice the gravitational softening length (|Springel et al.l 
bOOSi ). For radii greater than this convergence radius and 
less than the radius where the density of bound subhalo 
mass drops to ~ 80% of the local total mass density, the 
subhalos are well fit by an "Einasto" profile. 



Hp{r)/p-2] = {~2/a)[{r/r-,r-l]. 



(3) 



Here p_2 and r_2 are the scale density and scale radius (at 
the point where the density profile has the isothermal slo pe) 
respectively, and we take a = 0.17 l|Navarro et al.ll2010D . In 
our kinematic analysis, we extrapolate an Einasto fit to each 
subhalo when it is necessary to evaluate the mass distribu- 
tion at radii < 100 pc, and we use the directly determined 
mass profiles at all larger radii. 
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To model the three-dimensi onal stellar density profile, 
p-k, we use functions of the form (|Zhaoll 19971 ): 



p*(r) oc 



(4) 



where x = r/ro and {a,b,c,ro} are free parameters that 
will be estimated in the next section by fitting the observed 
surface density profile of each satellite. We focus on cuspy 
central profiles (0 ^ a ^ 1) because we find that they are 
required to fit the observed, nearly flat velocity dispersion 
profiles if we assume isotropic stellar velocity dispersions 
{P{r) = 0) and an Einasto halo profile. Although such cuspy 
profiles have not been used previously in studies of the MW 
satellites, they are, in fact, required to fit the inner surface 
brightness profiles of elliptical galaxies of all luminosities, 
including faint ones iGebhardt et al.lll996l ) and so seem a 
priori quite plausible for dSph galaxies also. 

If the stellar mass is everywhere negligible compared to 
the dark matter mass, the projected velocity dispersion in 
Eq. [2] is independent of the constant of proportionality in 
Eq. [4]that sets the stellar mass-to- light ratio, M*/L*. How- 
ever, if the stars contribute significantly to the potential of 
the galaxy, then we must determine the appropriate normal- 
izing factor for Eq. [4]and thus its contribution to the overall 
mass distribution of the galaxy. For each of the dSphs we 
will t ake M^^/L» = 1, consistent with t he observational re- 
sults l|Matedll99AlColeman et al.ll2005l '). We find that smaU 
variations in M*/L*, indicative, perhaps, of multiple stellar 
populations or differences in stellar initial mass function, 
have little effect on the results we present below. 

As also noted above, throughout our analysis we will 
assume locally isotropic velocity distributions, j3{r) = for 
all r. This is a strong assumption and it is thus remarkable 
that we find that we can fit all the kinematic data without 
relaxing it. 



2.2 Velocity Distributions 

The preceding discussion demonstrates the well-known fact 
that the observable quantities and af^siR) are insuf- 

ficient to determine the mass profile, M{r), of a spherical 
system unless the velocity anisotropy, /3(r), is specified. Ad- 
ditional kinematic information is contained in higher order 
moments of the line-of-sight velocity distribution, so appro- 
priate modeli ng of these moments may constrain /3(r) and 
so M{r) fe.g. lGerhardlll993l : lLokas fc Mamonll2003l ). A fully 
consistent dynamical model must clearly match the full line- 
of-sight velocity distribution at all radii. 

If we assume /3(r) — 0, it is possible to invert the observ- 
ables h{R) and o-f^siR^) to obtain not only a unique M(r) 
but also the unique distribution function, /(e), which repro- 
duces these observables within the potential corresponding 
to M(r). This distribution function then determines the full 
line-of-sight velocity distribution at each R. Thus, once we 
have found a subhalo with M{r) consistent with the I*{R) 
and o-foaiR) measurements for a particular dSph, we can 
check the consistency of the resulting model by comparing 
its line-of-sight velocity distributions with those observed. 

To obtain these velocity distributions, we begin with 
the Eddington inversion formula. 



(5) 



where e = *I'(r-) + «^/2 is the binding energy, is the grav- 
itational potential, and v is the modulus of the velocity. 
Potentials for the stars and the dark matter can be sep- 
arately constructed numerically via the Poisson equation, 
V^^i = iirGpi. The indices on potential and density repre- 
sent a specific component, the dark matter or the stars. The 
total potential is then the sum of the two. 

The Eddington formula in Eq.[5]determines the velocity 
distribution as a function of the binding energy. However, to 
compare to observations we need line-of-sight velocity dis- 
tributions for a set of circular annuli. Defining vios as the 
component of velocity along the line-of-sight and perform- 
ing the appropriate weighting over three-dimensional radii r, 
the distribution of line-of-sight velocities at projected radius 
R is given by 



f{vios;R) oc 



rdr 



i?2 



(6) 



9{r) + vf J2 



where rios is defined by 2'I'(r-;os) ~ wfoj,; for a given velocity 
vioa, we determine rioa via a numerical root-finding algo- 
rithm. The normalization of Eq.|6]will not be important for 
the purposes of our discussion. 

Once we have determined the velocity distribution in 
Eq.|6]it is straightforward to construct higher order moments 
of this distribution. In particular, the n*'' moment of the 
distribution is given by 



Wloa\R) 



J Vroaf{vios;R)dVloa 
J fiVloa; R)dvioa 



(7) 



where, as in Eq. (6] we have explicitly written / as a func- 
tion of the line-of-sight velocity. As an example that will 
be important for us below, the RMS velocity determined 
from Eq.His a/K~R)- We are thus able to check our nu- 
merical calculation of the velocity distribution function by 
comparing the RMS velocity determined from Eq. [7] to the 
equivalent quantity determined from Eq.[T] 

Eq. |6] gives the theoretical velocity distribution at R, 
but in practice, to compare to the observations, we must de- 
termine the distribution of / at the position of each observed 
star, and then average it over all the stars in each annulus. 
Thus, our mean fiyioa) for a given annulus is the mean of 
the values at the position of the stars, with the individual 
/ distributions all normalized to unity. With this procedure 
there are no approximations related to finite bin size when 
comparing our theoretical model to the data. 



3 DATA ANALYSIS 

In this section we perform our analysis of the photomet- 
ric and kinematic data. We first discuss how we use the 
star count data to identify parameter values for Eq. |4]that 
describe each satellite well. We then describe our handling 
and interpretation of the line-of-sight velocity data, and the 
way in which we use these data to identify specific Aquar- 
ius subhalos that could host each satellite. In particular, we 
describe the criterion by which we judge goodness-of-fit for 
a given satellite-subhalo match. 
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3.1 Photometry 

We use photometric data from the following sources for 
the surface density profiles of our five sa^laxies: F ornax 
llColeinan et a l.' 2005); Sculptor (Bat tagha et al.ll2008l); Sex- 
tans ( Irwin fc Hatzidimitriou .1995. 1: Leo I ( Smolcic et al.l 
|2007^ : Carina jMunoz et al-lbOOd ). Traditional fits to these 
data sets have used King or Plummer models, the latter 
corresponding to {a, h, c} = {0, 2, 5} in Eq. 2] Such fits typi- 
cally fai l to reproduce the measured s tar counts in the outer 
regions (jlrwin fc Hatz idimitrioullTggsl '). In addition, they en- 
force a constant density core which is consistent with star 
counts in some globular clusters but not with photometry of 
the inner regions of brighter early type galaxies, almost all 
of which show inner cusp s corresponding to a v alues signifi- 
cantly greater than zero iGebhardt et al.lll996l '). 

We take {a,6, c, ro} in Eq. |3]as free parameters to be 
adjusted when fitting the observed star count profiles. We 
perform a standard Abel projection of p*{r) to obtain I*{R) 
and we determine the free parameters for each satellite via 
a standard minimization procedure. In performing these 
fits, we find that there is a complex degeneracy in the space 
spanned by the four parameters. Motivated by reasons that 
we discuss in detail in Section[4l we focus on models in which 
the 3D stellar profiles are characterized by a central shallow 
cusp and a relatively sharp turnover to a steep outer power 
law. For Sculptor, Carina, and Sextans we will specifically 
adopt a central cusp with a — 0.5, and b = 3, while for 
Fornax, we will use a — 1 and 6 = 4. In Leo I, the star counts 
do appear to require a core, and we adopt a = and a similar 
transition to the steep outer power law. With a and 6 fixed 
a priori, we vary the remaining parameters, {c, ro}, in order 
to minimize x^- In ^11 cases this results in reduced values 
near unity, indicating an acceptable fit, and also within the 
90% c.l. of the minimum values attainable by varying all 
four parameters independently. The resulting parameter sets 
are given together with the corresponding /{Nbin — 5) in 
Table [T] Note that since our goal is to demonstrate that the 
observations are consistent with simple spherical, isotropic 
models within ACDM subhalos, it is not necessary for us to 
choose the best-fit profile parameters; rather we need only 
show that the parameters we do choose are consistent with 
the star count data. 

In Figure[T]we plot these surface density profiles for each 
satellite on top of the observed data. The scale radii vary 
over the range 0.29 kpc (Carina) ^ ro ^ 0.67 kpc (Fornax) 
and the outer slopes over the range 3.3 (Sextans) ^ c ^ 
7.5 (Leo). As noted above, the degeneracies allow signifi- 
cant variations in these quantities, particularly if a and 6 
are allowed to vary away from the values we have chosen. 
Our choices are motivated in part by simplicity (e.g. for a), 
in part by experimentation, determining which parameter 
ranges allow good fits also to the kinematic data (see be- 
low). For Sextans such considerations lead us to settle on 
a relatively shallow outer density profile, while for Leo the 
data force us to a steep outer profile. Note that in all cases, 
the star counts were actually carried out in elliptical annuli. 
The radial coordinate plotted is the geometric mean of the 
major and minor ajces which we expect to correspond best 
to the count profile for circular ann uli. (The typical elliptic- 
ities of these satellites are ~ 0.3 l|lrwin fc Hatzidimitrioul 
[l99i)). 
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Figure 1. Projected surface density profiles for each of the satel- 
lites that we consider, fit using the formula in Eq. (4] The values 
chosen for the parameters (a, b, c, ro) for each galaxy are given in 
Table [Tl together with the corresponding P^r degree of free- 
dom. Our procedures for selecting these parameters are outlined 
in section [3] where we also give references for the observational 
data plotted in each panel. 



3.2 Kinematics 

The kinematic datasets that we use co nsist of line- o f-sight 
stellar v elocities f r om t he samples of iMateo et al.l (|2008h 
and Wal ker et al] (|2009l). The latter use an "Expectation 
Maximization" method for evaluating membership and re- 
moving contaminants fr om each sample, an d we consider 
only those stars for which IWalker et al.l (|2009| ) assign > 90% 
probability of membership. The resulting numbers of stars 
are listed in Table [1] For Leo I, which is the only galaxy 
in our sample with out published mem bership probabilities, 
we use data from iMateo et al.l l|2008l ). and consider those 
stars as members that have velocities in the range from 
240 to 320 km s~^. As this range of velocities is well sepa- 
rated from that of MW foreground stars, it is unlikely that 
this sample suffers significant contamination. Other meth- 
ods for cleaning dSphs fro m contaminating MW ha lo stars 
have been considered (e.g. iKlimentowski et al.ll2007l l: these 
typically reduce the velocity dispersion at outer radii. The 
I Walker et all (|2009i l membership cuts appear appropriate 
for our analysis here. 

For each satellite, we bin the velocity data in a series 
of circular annuli and estimate the mean square line-of-sight 
velocity in each annulus as 



6-^ = (u^) - (e^). 

Here we define the velocity of a star as 



■ Vo,, 



(8) 

■Vo, where 
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Vo,i is the observed velocity of the j"* star and Vo is the mean 
of these velocities over all stars in the galaxy. The quantity 
represents the measurement uncertainty of the i**^ star, 
and angle brackets represent an average over all the stars in 
a radial bin. We further assume that the error on Vo is neg- 
ligible and that the actual velocities are uncorrelated with 
their measurement error. With these assumptions, is an 
unbiased estimator of the corresponding population quan- 
tity, and approximating the sampling distributions of («^) 
and (e'^) as normal, the uncertainty on a can be estimated 



2N {w2) - (e2) ■ 



(9) 



Given an estimate of the intrinsic velocity dispersion 
profile of each satellite based on Eq. [8l we step through all 
the subhalos in the six Aquarius simulations to determine 
which subhalo has the (spherically averaged) potential that 
best describes the data. Specifically, for each Aquarius sub- 
halo, we derive a spherical potential from the mass profile 
M(r) and then use the Jeans equation ([2]) to calculate the 
line-of-sight velocity dispersion profile, aiosiR), which cor- 
responds to the model star count profile of Table [T] and an 
everywhere isotropic velocity dispersion tensor. This line-of- 
sight velocity dispersion is then averaged over the positions 
of all the stars in each annulus to predict the population 
mean square velocity within that annulus. For each satellite- 
subhalo pair we then determine the quantity 



2 

X = 



E[(Ti — (Jlos{Rx) 

1=1 



(10) 



where A^bins is the number of annuli and is the mean value 
of the projected radius of the stars in the i**" annulus. For a 
given satellite, it then follows that the best fitting Aquarius 
subhalo is the one that minimizes Eq. 1101 

Once a "best" subhalo has been identified in this way, 
we can quantify whether it actually provides an acceptable 
fit by comparing the value from Eq.[TO]to the theoretical 
distribution of for A^'bins degrees of freedom. If p is the 
fraction of the theoretical distribution at larger values than 
the measured x^i then we can exclude the hypothesis that 
the observed satellite has isotropic velocity dispersions and 
is hosted by this "best" subhalo at confidence level 1 — p. 
(Note that, given our assumptions, there are no free parame- 
ters when comparing observed and predicted dispersion pro- 
files for a specific subhalo.) If p is not very small, then we 
conclude that the observed satellite could be hosted by a 
ACDM subhalo. Note that the converse does not apply. If 
p is very small, the observed satellite could still live in a 
ACDM subhalo if it has significant velocity anisotropics. 



4 RESULTS 

In this section we turn to the implementation of the algo- 
rithms described above. We begin by finding the Aquarius 
subhalo that best matches the line-of-sight velocity disper- 
sion of each satellite under the assumption of negligible ve- 
locity anisotropy and for the model stellar density profile we 
have fitted to the observed counts. We then check whether 
the line-of-sight velocity distributions of these models are 
consistent with those observed. 
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Figure 2. Line-of-sight velocity dispersion for our five satellites. 
The solid curves show the dispersion predicted by inserting the 
potential determined from the best fitting Aquarius subhalo and 
the photometric profile of Table [Tl into Eq. [T] assuming no veloc- 
ity aniso tropies. The sym bols show th e observational data taken 
from Ma teo et al] llooi) (Leo I) and IWalker et al] llooi) (For- 
nax, Carina, Sculptor, and Sextans). The errors on the velocity 
dispersion in each bin are assigned according to Eq. [9] 



Table 1. Number of member stars with measured radial velocities 
in each of our five galaxies, together with the parameters in Eq. |4] 
for our preferred fits to their star count profiles, as shown in Fig.[l] 
The final column gives the value of per degree of freedom for 
these count profile fits. 



Satellite 


# of stars 


a 


b 


c 


TO [kpc] 


xVd.o.f 


Fornax 


2409 


1 


4 


4.5 


0.67 


1.0 


Leo I 


328 





3 


7.5 


0.40 


1.6 


Carina 


758 


0.5 


3 


5.3 


0.29 


1.1 


Sculptor 


1392 


0.5 


3 


5.5 


0.32 


0.4 


Sextans 


424 


0.5 


3 


3.3 


0.44 


0.1 



4.1 Best-fitting subhalos 

Figure [5] compares the observed velocity dispersion profiles 
of our five satellites to those predicted by Eq. [2] when a 
stellar system with a star count profile given by Eq. |4]with 
the parameters in Table [T] with a stellar mass-to- light ratio 
of 1, and with negligible velocity anisotropy, is embedded in 
the Aquarius subhalo that fits best according to the criterion 
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Figure 3. Circular velocity profiles for the best-fitting subhalos 
(black squares) and the next-best-fitting subhalos (red diamonds) 
for our five satellites. The symbols give [GM{r)/r)]^^'^ which was 
used to calculate the potential from the dark matter component 
in the velocity dispersion profiles of Fig. [2] The profiles of each 
pair of subhalos are almost identical over the range encompassed 
by the stellar kinematic data. At larger radii deviations can be 
large. 



of Ea llOl The p- values for these best-fitting subhalos are 0.6, 
0.5, 0.6, 0.2 and 0.8 for Fornax, Leo I, Carina, Sculptor, and 
Sextans, respectively. Thus in all cases the fit appears good 
from a statistical point of view. 

It is important to note that when matching the data in 
Fig- El the photometric parameters chosen in Eq. U strongly 
affect the quality of the fit. For example, for Fornax we are 
able to find an Aquarius subhalo that matches both the 
photometry and the kinematics with an acceptable p-value 
only if a > 0.8. This motivates choosing a = 1 for this galaxy. 
For Sculptor, Carina, and Sextans the data constrain the 
central slope of the stellar density profile less strongly, so for 
these satellites we choose a = 0.5. As noted above, the star 
count data for Leo I force a value of a close to zero. Hence 
we pick a — for simplicity, although a slightly negative a 
gives a somewhat better fit to the velocity dispersion data. 

In Figure [3l we show circular velocity curves for the two 
"best" subhalos for each of the satellites; the black squares 
correspond to the subhalos plotted in Fig. [21 There is es- 
sentially no difference between these pairs of curves at radii 
where we have kinematic data. This refiects the fact that. 



for our assumptions, the potential can be derived directly 
from the observational data so any "acceptably fitting" sub- 
halo will have to resemble the result of this exercise quite 
closely. At larger radii the potential is effectively uncon- 
strained, however, and the profiles of the two subhalos can 
differ dramatically. For Fornax, Leo I, Carina, Sculptor and 
Sextans the maximum circular velocities of the best-fitting 
subhalos are 21, 22, 16, 26, 12 km s"\ respectively. The 
corresponding present-day dark matter masses are 7, 2, 2, 
15, and 1 x 10* M©. 

These results allow us to conclude that the Aquarius 
simulations include at least one subhalo which is an accept- 
able host for each of the five satellite galaxies, even under 
the restrictive assumption of negligible velocity anisotropy. 
We can extend this analysis to estimate the total number of 
Aquarius subhalos that are statistically consistent with the 
kinematic and photometric data for each satellite. Specifi- 
cally we count those subhalos that have a p-value exceed- 
ing 10%. (Remember we are still enforcing isotropic velocity 
distributions and the specific profile parameters of Table [T]). 
For Fornax we find 13 acceptable subhalos among the six 
simulations; for Sculptor there are 16 such subhalos; for Ca- 
rina there are 8 subhalos; for Leo I there are 37; and for 
Sextans we find over 100 acceptable hosts. For Sextans this 
large number of hosts is a reflection both of the fact that its 
velocity dispersion is lower and that the errors on the data 
are relatively large. 

4.2 Full velocity distributions 

The preceding analysis shows that ACDM subhalo poten- 
tials can fit the star count and line-of-sight velocity disper- 
sion profiles of Milky Way satellites to good accuracy, even 
under the assumption of negligible velocity anisotropy. We 
now use the techniques described in Section [2] to study if 
these same models provide good descriptions of the full line- 
of-sight velocity distributions, for example, of their higher 
order moments. 

The line-of-sight velocity distribution at projected ra- 
dius i? is a convolution of the intrinsic distribution from 
Eq. [6] with the measurement uncertainty of the stellar ve- 
locities. The former has to be averaged over the radial po- 
sitions of all the stars in each of our annuli, and we model 
the latter by a Gaussian with dispersion equal to the mean 
quoted velocity error of these same stars. This measurement 
uncertainty ranges from ~ 3 km s^^ for all the Carina bins 
and the innermost Sculptor bins, down to ~ 1.5 km s^^ for 
all the Fornax bins. 

The resulting line-of-sight velocity distributions are 
shown for each of our five galaxies in Fig. [l] For each galaxy 
we used four circular annuli containing nearly equal num- 
bers of stars. Our motivation for this binning scheme is to 
retain some information about the variation of distribution 
shape with radius while keeping a large number of stars in 
each bin in order to better constrain the shape of the dis- 
tribution. This results in fewer bins in Fig. [4] than in Fig. [2] 
for Fornax, Sculptor, and Carina. In each panel of Fig. 2] 
the solid curve is the line-of-sight velocity distribution cal- 
culated as above and normalized to have the same area as 
the corresponding histogram. 

To quantify the level of agreement between model and 
data in Fig. |4l we compare the distributions of \vi\ in each 
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Figure 4. Velocity distributions in four bins for eacli of our satellites. In each panel, a solid curve shows the theoretical distribution 
averaged over the radial positions of all the stars in the bin, and then smoothed with a Gaussian representing the typical observational 
error on the stellar velocities. Labels at the upper right list both the total number of stars in each bin and the approximate radial range 
they encompass. 



panel using a KS test. We take the modulus here because 
the distribution of the line-of-sight velocity relative to the 
galaxy mean is expected to be symmetric about zero for 
any equilibrium model (even rotating and/or non-spherical 
ones) after averaging over a circular annulus. As a result, 
all shape information is contained in the distribution of \vi\, 
and restricting the test in this way enhances its sensitivity 
to higher order moments. The maximum difference between 
the normalized cumulative distributions of \vi\ for data and 
model is then a measure of the confidence level at which 
we can reject the null hypothesis that our simple isotropic, 
spherical model represents the full, observed line-of-sight ve- 
locity distribution in the annulus. 

Results of this KS test for each of the four annuli and 
for all of our satellites are shown in Table [21 with bins 1-4 
ordered by increasing radius. These values indicate that our 
predicted velocity distributions are generally in good agree- 
ment with the data. The annuli with the lowest probabilities 
are bins 3 and 4 of Fornax and bin 4 of Leo I; for Fornax 
bin 3 and Leo I bin 4 the null hypothesis can be excluded 
with > 99% confidence. As can be seen in Fig. IH the mea- 
sured line profiles appear less peaked (platykurtic) than the 
models in these annuli. For Fornax the effect is quite weak, 



but is nevertheless significantly detected because of the large 
number of stars involved. For all the other panels of Fig. [4] 
differences are less significant and are quite small. Note that 
the kurtosis of these distributions is expected to be quite 
sensitive to velocity anisotropy, so the fact that our models 
fit fairly well can be taken as an indication that anisotropies 
are probably weak. 

A more direct measurement of the kurtosis of the line- 
of-sight velocity distributions can be obtained by estimating 
their fourth moment directly. For each annulus, we calculate 
a sample kurtosis from the A'^ stars it contains as 

^ = {v')/{v'f-3. (11) 

The kurtosis is defined so that a Gaussian model gives k = 0. 
W e appr oximate the uncertainty of the sample kurtosis by 
y^24:/N, the scatter expected for random samples from a 
normal distribution. For our theoretical model, we calculate 
second and fourth moments from Eq. [71 after smoothing 
to account for measurement errors, and we then substitute 
these into Eq. [TTlto obtain the predicted kurtosis. 

The results of this exercise are shown in Fig. \S\ The 
models predict very little kurtosis in almost all annuli, and 
the data agree with this in most cases. Comparing Fig. [5lto 
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Table 2. KS probabilities for the maximum difference between 
the observed and modeled cumulative distributions of \vi\ within 
four equally populated annuli in each of our observed satellites. 
Bins 1-4 correspond to the annuli of Fig.[4]ordered from inside to 
outside. The complement of each of these values represents the 
confidence level at which the hypothesis that the data are drawn 
from the theoretical distribution can be rejected. 



Satellite 


bin 1 


bin 2 


bin 3 


bin 4 


Fornax 


0.05 


0.15 


0.003 


0.02 


Leo I 


0.46 


0.52 


0.31 


0.003 


Carina 


0.51 


0.44 


0.29 


0.34 


Sculptor 


0.32 


0.68 


0.62 


0.67 


Sextans 


0.41 


0.81 


0.03 


0.97 



the second moment data in Fig. [2] clarifies the origin of the 
discrepancies uncovered by our KS test. For example, in the 
outer annuli of Fornax, the model distribution is (slightly) 
narrower than the data while the kurtosis estimates are al- 
most in agreement. For the 4th annulus of Leo I, the model 
RMS is again smaller than for the data, and the data are 
more platykurtic. For Leo I, the 2nd annulus also shows a 
marginally significant platykurtic signal, while for Sextans 
bin 3 the data are more leptokurtic than the model. How- 
ever, in the latter two instances the number of stars is too 
small for the KS test to indicate a significant discrepancy. 
The rather large bin-to-bin fluctuations in the observational 
estimates of kurtosis suggest that our Gaussian error bars 
may be underestimating the true sampling uncertainties. 

From this analysis, we conclude that the observed line- 
of-sight velocity distributions agree surprisingly well with 
our spherical, isotropic models embedded in ACDM sub- 
halo potentials. The remaining differences can plausibly be 
ascribed to departures from isotropy, from spherical sym- 
metry, from dynamical equilibrium, or (more likely) from a 
combination of these. Detailed modeling is, of course, nec- 
essary to test this possibility, but it is beyond the scope of 
the present paper. Our main goal here has been to show 
that simulated ACDM halos contain subhalos with poten- 
tials consistent with those hosting the observed stellar pop- 
ulations of Milky Way satellites, even if these are assumed 
to have negligible velocity anisotropy. 



5 DISCUSSION AND CONCLUSION 

We have investigated whether the gravitational potentials 
of subhalos in N-body simulations of ACDM halo forma- 
tion are consistent with the high-quality photometric and 
kinematic data available for five of the brighter satellites of 
the Milky Way. We find that a direct mapping is, in fact, 
possible between each of these satellites and a subset of the 
dark matter subhalos in the six high resolution simulations 
of the Aquarius Project. Star count profiles with inner cusps 
scaling as r~" with ^ a ^ 1 can provide good fits to the 
observed counts. Placed in the measured Einasto-like po- 
tentials of appropriately selected subhalos, they also fit the 
observed, nearly flat line-of-sight velocity dispersion proflles 
very well, even under the restrictive assumption of negligi- 
ble velocity anisotropy. Such isotropic models fit the shapes 
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Figure 5. Comparison of the predicted and observed kurtosis of 
the line-of-sight velocity distribution in four equally populated 
annuli in each of our five satellites. The kurtosis is defined to be 
zero for a normal distribution. The black squares are direct obser- 
vational estimates for each annulus, with errors given by 24/7V, 
where A'^ is the number of stars. The red diamonds are theoreti- 
cal predictions from our isotropic models after convolution with 
the measurement errors and averaging over the stellar positions 
in each annulus. 



of the observed line-of-sight velocity distributions well, in 
addition to their second moments. 

We have measured the present-day maximum circular 
velocities of the "best-fit" subhalos for each of these five 
satellites. These range from 10 to 30 km s~^. Subhalos con- 
sistent with hosting the observed systems at 2 = have 
peak circular velocities (i.e. the largest maximum circular 
velocity they ever had) ranging from 12 to 50 km s"\ The 
maximum past masses of the main progenitors of these sub- 
halos range up to ~ 5 x 10^ Mq , while their masses at 2: = 7 
(the ap proximate lower b ound on the redshift of reioniza- 
tion (jbunklev et al.ll2009l ')') range up to ~ 10** M©. At z = 
their Galactocentric distances range from 40 to 400 kpc. 

Our results indicate that current data on faint Milky 
Way satellites are consistent with these galaxies living in 
ACDM halos. They do not, however, explain why galaxies 
living in such subhalos should have the observed properties. 
Expl oring this issue is an import ant task for future work 
(e.g. iLi et aLlbOlOl: ICooper et al.r 2009: Sa wala et al.ll2010l : 



lOkamoto fc Frenkl(2009r 



Busha et al.ll201Gl ). 
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